Competing itinerant and localized states in strongly correlated BaVSa 
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The electronic structure of the quasi- lowdimensional vanadium sulfide BaVSs is investigated for 
the different phases above the magnetic ordering temperature. By means of density functional theory 
and its combination with dynamical-mean field theory, we follow the evolution of the relevant low- 
energy electronic states on cooling. Hence we go in the metallic regime from the room temperature 
hexagonal phase to the orthorhombic phase after the first structural transition, and close with 
the monoclinic insulating phase below the metal-insulator transition. Due to the low symmetry 
and expected intersite correlations, the latter phase is treated within cellular dynamical mean- 
field theory. It is generally discussed how the intriguing interplay between band-structure and 
strong-correlation effects leads to the stabilization of the various electronic phases with decreasing 
temperature. 
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I. INTRODUCTION 

Since its first characterization^ in 1969, the under- 
standing of the complex electronic structure of the vana- 
dium sulfide BaVSs poses a longstanding problem in con- 
densed matter physics^i^i^i^i^. Numerous experimental 
and theoretical studies have revealed a delicate coupling 
between orbital, spin and lattice degrees of freedom over 
a wide temperature range. On cooling, BaVSs exhibits 
three continuous phase transitions, starting with a struc- 
tural hexagonal-to-orthorhombic transition at Ts~240 K 
in the metallic regime. The latter vanishes at ~70 K 
where a metal-to- insulator transition (MIT), accompa- 
nied by a lattice transformation from orthorhombic to 
monoclinio^i^, to a still paramagnetic phase takes place. 
A final magnetic transition marking the onset of an in- 
commensurable antiferromagnetic order— takes place 
at Tx- 30 K. 

The underlying driving forces for these transitions and 
the specific nature of the respective phases is to a large 
extent still a matter of debate. It was showr^ that the 
MIT may be driven to zero temperature at high pressure, 
and the supression of the insulating phase leads to non- 
Fermi-liquid and quantum-critical behavior^^. This ob- 
servation is not only adding even more complexity to the 
already existing problems, but also underlines the tricky 
nature of the electronic structure. At ambient pressure 
the MIT is early announced by strong precursive behav- 
ior such as a large increase of the Hall coefficient^ and 
a wide one-dimensional (ID) lattice-fluctuation regimei^ 
along the c axis of the system. In fact, it seems to be es- 
tablished^i^ii^ii^ that the MIT may be described in terms 
of a charge-density wave (CDW) instability. However, 
BaVSs is not a textbook Peierls system. The dc con- 
duction anisotropy is rather smalU^ (<Jc/cra~3-4), and 
the "metallic" phase above Tmit displays a high resistiv- 
ity (a few ml^cm) and metallic-like behavior {dp/dT>0) 
only above a weak minimum at ~150 K, below which 
it increases upon further cooling^ii^. Moreover, local- 
moment behavior is revealed from the magnetic suscep- 



tibility, with an effective moment of approximately one 
localized spin-1/2 per two V sites. At Tmit, the sus- 
ceptibility rapidly drops, and the electronic entropy is 
strongly suppressed^^. 

Because of the nominal V^+ valence, BaVSs belongs to 
the family of 3d^ compounds with a t2g manifold span- 
ning the low-energy sector. In the hexagonal phase the 
latter consists (per V ion) of an Aig and two degen- 
erate Eg states. The remaining Cg states are strongly 
hybridized with the S(3p) states and have major high- 
energy weight. Below Ts the degeneracy of the EgS is 
lifted in the orthorhombic phase. In both phases the 
primitive cell includes two formula units, where the V 
and S ions form chains of face-sharing VSs octahedra 
along the c axis. The intrachain V-V distance is less 
than half the interchain distance. It follows^ that the 
Aig orbital is mainly directed along the chain, forming 
a broader band due to the significant overlap of neigh- 
boring intrachain V ions. On the contrary, the lobes of 
the Eg orbitals point inbetween the sulfur ions, i.e., do 
not hybridize strongly with their environment, leading 
to comparably narrow bands. These simple characteriza- 
tions hold essentially also for the monoclinic insulating 
phase, yet the primitive cell is doubled and the resulting 
four V ions in the basis are now all inequivalent by sym- 
metrys^iii. The CDW mechanism has led to a tetramer- 
ization, yet no evident charge disproportionation among 
them was detected^^. The measured charge gap^ii^ii^ii^ 
of about 40 meV is twice as large as the apparent spin 
gap2U, pointing once more towards the relevance of elec- 
tronic correlations^!. 

In Ref. p^l22| the orthorhombic phase above the 
MIT was investigated, and it was argued that strong 
electronic correlations are responsible for a substantial 
charge transfer within the t2g states, leading also to im- 
portant Fermi-surface changes in comparison to a weak- 
correlation treatment. Here we go further by tracing the 
low-energy states of BaVSs all the way from room tem- 
perature down to 40 K (just above the final magnetic 
transition). Although the local environment of the V site 
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does not change dramatically, the electronic structure ap- 
pears to be rather sensitive to the temperature changes. 
This originates from the subtle balance of kinetic energy 
versus Coulomb interaction in the electronic system, a 
characteristic of strongly correlated materials. In fact, 
this vanadium sulfide presents an interesting realistic re- 
alization of one of the basic problems in strongly corre- 
lated physics: there is nominally one electron in the low- 
energy sector and two distinct orbital states, one forming 
a broader and two forming narrower bands. Hence de- 
pending on temperature, nature shall find the best com- 
promise between kinetic energy gain and potential cost 
due to mutual Coulomb interaction in this multiorbital 
scenario. 



II. THEORETICAL FRAMEWORK 

For the investigation of competing band-structure 
and many-body effects in realistic materials, the com- 
bination^^'^^ of density functional theory (DFT) and 
dynamical-mean field theory (DMFT) has recently 
proven to be a powerful approach. 

For the DFT part we used the local density approxi- 
mation (LDA) to the exchange-correlation energy. The 
corresponding calculations were performed with a mixed- 
basis pseudopotential code^^. It uses normconserving 
pseudopotentials, and plane waves supplemented with 
some few non-overlapping localized functions in order to 
represent the pseudo crystal wavefunction. 

Since the low-energy physics of BaVSs is dominated 
by the t2g states, the so-called LDA+DMFT calculations 
were performed for the corresponding three-band subset. 
The latter was derived from the full band structure via 
the maximally-localized Wannier function (MLWF) con- 
struct ion^^i^i. Hence the local orbitals which form the 
impurity in the DMFT context stem from the associ- 
ated Wannier functions (WFs) and the low-energy LDA 
hamiltonian H{k.) is expressed with respect to these or- 
bitalsi^. To be specific, by making reference to the for- 
malism outlined in Ref. Hi, in all LDA+DMFT calcula- 
tions presented here the set of correlated orbtials C was 
identical to the set W of the WFs forming the minimal 
LDA hamiltonian. Hence the impurity Green's function 
is computed in DMFT for finite inverse temperature P 
via 

k 

(1) 

where ujn={'^n + are the Matsubara frequencies 

and 5] is the self-energy matrix for the strongly corre- 
lated orbitals. For the local interacting hamiltonian Hi^t 
the following representation restricted to density-density 



terms only was used: 

Hint = U^hm^hmi^ — ^ flmaflm'a 

m mm' a 

m^m' 

H" 2 / ^ '^ma'^m' a • \'^) 

mm' a 
m^m' 

Here fima—^Xna^ma^ whcrc m,cr denote orbital and spin 
index. The following parametrization of U' and U" has 
been proven to be reliabl o^^i^^ in the case of t2^-based 
systems: U'=U-2J and U"=U-ZJ. We utilized the 
quantum-Monte Carlo (QMC) formalism after Hirsch- 
Fye^^ to solve the impurity problem. 

For hexagonal and orthorhombic BaVSa there are only 
symmetry-equivalent V ions in the primitive cell and 
moreover interatomic correlation effects are not expected 
to be of crucial importance. Hence a single-site DMFT 
approach to describe the strong-correlation effects was 
employed for those phases. Thereby the inverse tempera- 
ture was always /3=30 eV~^ and the number of time slices 
equaled 128 for the QMC method. However, since the 
paramagnetic insulating regime of the low-temperature 
monoclinic phase is associated with a CDW state, this 
approximation appears inadequate. We thus used a clus- 
ter formalism (for recent reviews see e.g. Ref. (3lll32ll33| ) , 
namely the cellular DMFT (CDMFT) approach in a re- 
alistic context— More explicitly, the linear cluster 
formed by the four symmetry-inequivalent V ions, each 
one decorated with a three-orbital t2g multiplet, was 
identified as the impurity for the CDMFT scheme. This 
amounts to a self-energy matrix YiiiuOn) that is not only 
off-diagonal in the orbital indices, but also in the site in- 
dices within the cluster. Note however that inter cluster 
components of the self-energy are neglected. The lat- 
ter fact results in the breaking of translational symme- 
try when computing pair correlations for the end sites 
of our linear cluster. However such a cluster approach 
should still be sufficient to describe the major qualita- 
tive changes in the correlated electronic structure origi- 
nating from the CDW instability. Since the undertaken 
cluster investigation is numerically very expensive within 
QMC (to our knowledge one of the largest up to now per- 
formed in the framework of realistic cluster-DMFT) we 
chose /3=25 eV~^ and used 90 time slices. 

III. RESULTS 

Stoichiometric BaVSa transforms on cooling succes- 
sively to crystal systems with lower symmetry, giving 
a hint to the generally low ordering energy. Table [J 
summarizes the basic crystal data for the three different 
phases that are studied in this work. Both, the hexag- 
onal and orthorhombic phases are associated with the 
metallic regime, while the monoclinic phase corresonds 
to the insulating system. Here we only investigated the 
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FIG. 1: (color online) BaVSs in the Cmc2i structure. The V 
ions are shown as smaller (red/gray) spheres, the Ba ions as 
larger (blue/dark) spheres. 



TABLE I: Experimental crystal data used for the investiga- 
tion of BaVSs. 





RT 


100 K 


40 K 


crystal system 


hexagonal 


orthorhombic 


monoclinic 


space group 


PGs/mmc 


Cmc2i 


Im 


a (a.u.) 


12.71 


12.77 


12.78 


h (a.u.) 


22.01 (= aVS) 


21.71 


21.65 


c (a.u.) 


10.63 


10.58 


21.15 








90.045 


experiment 


Ref. 36 


Ref. 36 


Ref. 8 



paramagnetic phase of the insulator and excluded the 
magnetically ordered phase below Tx- 



A. LDA study of the metallic regime 

At room temperature (RT) BaVSs crystallizes in the 
hexagonal {P6s/mmc) structure^ with two formula units 
in the primitive cell. All symmetry operations of the 
hexagonal group apply to this structure, and all Ba, V 
and S ions in the cell form one single symmetry class, 
respectively. The V ions within the chains are aligned 
straightly. As noted by Mattheis^l, the variable param- 
eter x(S) which determines the inplane S-S distances 
deviates slightly from the ideal hexagonal value, i.e., 
x(S)=0.1656 while Xideai=l/6. This means that the in- 
plane equilateral S-S-S triangle above and below an V 
ion is contracted and the intrachain S-S bond lengths are 



decreased in comparison to the interchain lengths. 

At Ts the crystal system transforms from hexagonal to 
orthorhombic, leading to a structure with (Cmc2i) space 
group. The V chains are now zigzag distorted in the be 
plane of the lattice. Although still two formula units form 
the primitive cell, the symmetry class of the sulfur ions 
has split into two: both SI ions are positioned at (4a) 
apical sites on the b axis, while the four S2 ions occupy 
(8b) sites (see Fig. [1]). From Tab. [J it is seen that the 
structural transformation results in only minor changes 
in the lattice parameters, with a maximum 1.5% per cent 
contraction of the b axis. There are several LDA studies 
for BaVSs above the MIT^^^^. In addition to those, 
a direct comparison of the LDA low-energy electronic 
structure for the hexagonal and orthorhombic phases via 
Wannier construction for the t2g manifold is presented 
here. 

Figure [2] displays the LDA density of states (DOS) for 
the two crystal structures. In each case the dominance 
of the t2g states at low energy, with a prominent peak 
right at the Fermi energy, is evident. For the Aig band 
the dispersion is indeed reminiscent of ID characteristics, 
however below the Fermi energy features are changed due 
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FIG. 2: (color online) Total and local LDA DOS for BaVSs 
in (a) the PQs/mmc structure and (b) the Cmc2i structure. 
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FIG. 3: (color online) (top) LDA band structure and (bot- 
tom) derived t2g Wannier bands for PGs/mmc-BaVSa. Color 
coding for (a): Aig (blue), Egi/Eg2 (red). 



to hybridization with the S(3p) states. The weight of the 
latter states is reduced in the energy range [—1,0] for 
the Cmc2i structure compared to the P6s/mmc struc- 
ture. One might interpret this as some decoupling of Aig 
and S{3p) throughout the structural transition. The very 
large DOS at the Fermi energy, common to both struc- 
tures, due to the Eg states render instabilities towards 



broken-symmetry phases very likely. 

The LDA band structures in Fig. [3^ and Fig. 2^, with 
"fatbands" exhibiting the weight of the respective t2g 
states on the different bands^^, show the folded character 
of the Aig band due to the two-formula unit primitive 
cell. The Aig bandwidth is dominated by the disper- 
sion along r-(A,Z), i.e., along the c* axis. While for the 
higher-symmetry P6s/mmc structure the folding prop- 
agates gapless through A, the corresponding bands are 
separated at Z for the Cmc2i structure. This hybridiza- 
tion between Aig and Eg is obvious in the upper triangle 
A-H-L and Z-E-T. Clearly seen is the hybridization of 
Aig with S(3p) resulting in "jumps" of the Aig charac- 
ter between different bands in the lower triangle F-K-M 
and F-C-Y. For both structures the Aig band cuts the 
Fermi level close to the zone boundary, leading to a 
nearly filled lower part of the folded band complex. It 
was noted in former works^^*^ that for Cmc2i this filling 
renders a CDW instability within the Aig band along F-Z 
impossible. 

The Eg states form very narrow, in some regions even 
nearly dispersionless, bands right at ep- An important 
difference occurs between the Fermi surface (FS) of the 
two metallic phases (see Fig. [5j). Although in both cases 
the FS consists of two sheets, the one for P6s/mmc is 
entirely located at the zone boundaries with a dominant 
Aig sheet and smaller Eg pockets around L. On the con- 
trary the FS for Cmc2i shows as the first sheet a substan- 
tial Eg2 electron pocket centered at F and Egi pillar-like 
structures on the h* axis extending along c*. The FS is 
completed by the Aig sheet, now extending deeper into 
the Brillouin zone (BZ). Despite the latter observation. 




FIG. 4: (color online) (top) LDA band structure and (bottom) 
derived t2g Wannier bands for Cmc2i-BaVS3. Color coding 
for (a): Aig (blue), Egi (red) and Eg2 (green). Fig. [3] 
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FIG. 5: (color online) (a) hexagonal and orthorhombic Bril- 
louin zones as well as LDA Fermi surface for BaVSs in (b) 
the PGs/mmc structure and (c) the Cmc2i structure from 
different perspectives. 
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FIG. 6: (color online) t2g Wannier functions for BaVSs in the 
crystal- field basis of (a) the PGs/mmc structure and (b) the 
Cmc2i structure (b). 



this quasi- ID sheet is neither strongly flattened nor have 
both parts the proper distance for nesting with the ex- 
perimental CDW vector (q^''''P^=0.5c*). 

The dispersions according to the derived three-band 
hamiltonian on the basis of the maximally-localized pro- 
cedure^^iSi are shown in Fig. 1314b . Because of the en- 
tanglement of the V{3d) bands with the S{3p) ones, the 
dominant Ai^-like band does not coincide with the true 
LDA bands in this minimal model^^. Still such a three- 
band approach should carry the essential physics in the 
low-energy regime. Figure [6] pictures the corresponding 
V(t2^) Wannier orbitals in the crystal- field basis, i.e., 
with vanishing onsite hybridization. For both phases it 
is seen that the Egi orbitals leak out on the S2 ions, 
while the Eg2 orbitals have weight on the apical SI atoms 
(recall that (S1,S2) are symmetrically equivalent only in 
the hexagonal phase). The Aig orbital hybridizes with 
the 3p orbitals on both sulfur-ion types. That this hy- 
bridization is indeed weakened in the Cmc2i structure 
may be derived from the reduced spread of the Aig WF 
shown in Tab. [Ill Interestingly, the spread is now even 
below the values for the Eg states. Note also that the 



TABLE II: Wannier centers Hw and spread (r^) of the t2p-like 
MLWFs constructed from a (6x6x6) /c-point mesh. The po- 
sitions of the symmetrically equivalent V sites in cartesian co- 
ordinates (in a.u.) read for PQs/mmc: Rvi = (0.00,0.00,0.00) 
and Rv2= (0.00,0.00,5.30) ; as well as for Cmc2i: Rvi = 
(0.00,0.46,-0.01) and Rv2= (0.00,-0.46,5.28). The V(2) site is 
symmetry-related to the V(l) site by the symmetry operation 
C2^'*Rvi+0.5. In the following only the data for VI is shown. 



structure 


WF 


Kw - Rvi (a.u.) 


(r^) (a.u.^) 




Aig 


0.00, 0.00, 0.00 


18.62 


PQs/mmc 


Egl 


0.00, 0.00, 0.00 


17.10 




Eg2 


0.00, 0.00, 0.00 


17.10 




Aig 


0.00, 0.30, -0.16 


16.60 


Cmc2i 


Egl 


0.00, 0.19, 0.35 


17.55 




Eg2 


0.00, 0.56, -0.31 


17.53 



TABLE III: Hopping integrals between the t2g Wannier or- 
bitals of BaVSs in the crystal- field basis. The first value cor- 
responds to the PQ^/mmc structure and the second to the 
Cmc2i structure, respectively. The term '00 1' shall denote 
the hopping to the nearest-neighbor V site within the unit 
cell. One of the nearest-neighbor V ions in the ab plane is 
located at TOO', while '110' and '110' are closest V ions along 
a and 6, respectively. Energies in meV. 









Egl 


Egl 


Eg2 


Eg2 


Aig 


Egl 


Aig 


Eg2 


Egl 


Eg2 


000 


395 


423 


200 


210 


200 


236 




















00 1 


-587 


-511 


90 


44 


-90 


-12 





-146 














001 


-61 


-86 


5 


14 


5 


3 





7 














100 


-49 


-35 


8 


14 


-40 


-26 


18 


-7 


-32 


-14 


41 


14 


110 


-49 


-26 


-63 


-76 


31 


29 


-37 


-28 





-2 





-12 


110 


3 


1 





2 


-5 


-6 





-2 















WF centers are identical with the V positions only for 
the hexagonal structure, while there are some shifts for 
the orthorhombic structure. The main qualitative differ- 
ence in the hopping integrals (Tab. IIIip between the two 
structures is the emerging substantial hybridization be- 
tween Aig and Egi for Cmc2i. Note that generally the 
Eg hoppings are rather isotropic, whereas the dominat- 
ing hopping along the c axis in the case of Aig is obvious. 



B. LDA-hDMFT study of the metallic regime 

The LDA approach presented in the last section under- 
lies the assumption that the mutual interactions between 
the electrons may be cast into a static local exchange- 
correlation potential within an effective single-particle 
description. However it is known that for strongly cor- 
related systems such a description is likely to fail. Since 
one indeed expects rather strong correlations within the 
quarter-filled V(t2^) states of BaVSs, also due to the 
specific characteristics of very narrow Eg bands and the 
broader Aig band, we employed the LDA+DMFT frame- 
work to explicitly include many-body effects in the elec- 
tronic structure. 
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By identifying the derived t2g WFs in the crystal-field 
basis as the subspace of correlated orbitals, we derived 
the k-integrated (local) spectral functions p{uj) shown in 
Fig. [71 Clearly seen is the transfer of spectral weight 
from the quasiparticle (QP) peaks to lower/higher en- 
ergies in comparision to the local LDA DOS, especially 
for the Eg states. This corresponds to an inclusion of the 
atomic-like excitations important for states with substan- 
tial localized character, which is missing in the standard 
LDA picture. By varying T one observes^^ additionally 
that the strength of the Eg QP peak changes also more 
significantly. Thus the corresponding electrons are ef- 
fectively localized for a wide temperature range due to 
incoherence effects. 

Since the correlation effects infiuence the subtle ener- 
getic balance in this system, they are moreover responsi- 
ble for a substantial charge transfer between the relevant 
orbital a^^i^^ , resulting in different orbital-resolved fillings 
compared to the LDA ones (see Tab. |lVl). The LDA fill- 
ing of the Aig band is close to 70% in the hexagonal 
structure and does not change much for temperatures 
where the DMFT study was elaborated. A reasonable 
choice^^ for the Hubbard parameters U and J leads to 
a significant transfer of charge from Aig to Eg^ in order 
the overcome the large potential energy cost for occu- 
pying mainly the former orbital. The new fillings for 
Egi and Eg2 differ slightly due to a marginal hybridiza- 



TABLE IV: Orbital-resolved fillings for BaVSs from 
LDA+DMFT within the crystal-field Wannier basis. The 
QMC solver was used for T=390 K. 




FIG. 7: (color online) (right) Local spectral functions from 
LDA+DMFT in comparison to (left) the local LDA DOS for 
the t2g WFs in the crystal-field basis, for (a) the PQs/mmc 
structure and (b) the Cmc2i structure. The QMC solver was 
used for T=390 K. 



structure 


U,J (eV) 




Egl 


Eg2 


P63 / mmc 


0.0, 0.0 
3.5, 0.7 


0.67 
0.45 


0.16 
0.29 


0.16 
0.27 


Cmc2i 


0.0, 0.0 


0.59 


0.31 


0.10 


3.5, 0.7 


0.41 


0.45 


0.14 



tion between Aig and Egi in the Wannier hamiltonian 
already for the hexagonal structure. Note that Riling in 
LDA+DMFT is not equivalent to pure band filling in the 
LDA sense, because the atomic-like excitations are now 
also included. The essential change for the orthorhombic 
structure below Ts (which was here treated at same T 
within the QMC solver of DMFT) is the effective reduc- 
tion of the three-band to a dominant two-band problem. 
Due to the now substantial Aig-Egi hybridization the 
charge transfer is dominantly taking place between those 
two orbitals. Yet the overall occupation of the Aig or- 
bital is only little smaller than in the hexagonal phase. 
The Egl filling is now close to 50%, in good agreement 
with the experimentally observed local magnetic moment 
of about one free spin every other V ion. 

Besides the orbital-resolved filling, the change of the 
respective FS sheets in the metallic regime is also of large 
interest. However keep in mind that there is no straight- 
forward relation between those two issues, as only the 
total Fermi-surface volume is invariant when turning on 
correlations and no unique rule of how the individual 
sheets have to change can be derived. But this change 
can of course be calculated, which was done for the or- 
thorhombic phase with the Wannier basis in Ref. 14, and 
the change of the FS sheets is indeed in line with what 
one expects from the overall charge transfers (s. Fig. [8]). 
More explicitly, a strong Fermi-surface deformation was 
revealed, placing parts of the renormalized Aig sheet 






FIG. 8: (color online) comparison between (a) the LDA FS 
of orthorhombic BaVSs and (b) the corresponding QP FS 
derived from LDA+DMFT. 
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now in reasonably good position for a possible nesting 
with the experimentally determined q vector. An impor- 
tant observation was that the nesting should mainly take 
place away from the high-symmetry directions in the BZ, 
i.e., the Aig band along the F-Z direction should not be 
strongly involved in the direct nesting. As pointed out, 
one expects the Eg states to be essentially localized for 
higher temperatures, rendering the definition of a sharp 
FS rather difficult. For this reason, no renormalized FS 
was computed for the hexagonal phase. Due to the miss- 
ing Aig-Egi hybridization one would however expect that 
for this phase the overall Aig sheet is shifted more or less 
coherently in the BZ, contrary to the orthorhombic case. 



C. LDA study of the paramagnetic insulating 
regime 

BaVSs below Tmit is insulating with a monoclinic Im 
structure involving four formula units in the primitive 
cell^. The system is described as a CDW state with a 
dominant 2k f distortion^. Figure [9] displays the four in- 
equivalent vanadium ions along the chain together with 
numbers indicating the shift of the atomic positions with 
reference to the Cmc2i structure above the MIT (at 
T=100 K). It is seen that the shifts for V along the chain 
are relatively small, on the scale of 4% at most. Nonethe- 
less, a dominant disortion pattern for this CDW state 
may be identified. The mainly shifted ions are VI and 
V3 in our notation, whereas V2 and V4 only marginally 
change their positions. Since VI and V3 are shifted to- 
wards each other, with V4 inbetween, the tetrameriza- 
tion appears as an effective trimerization, isolating the 
V2 site. As a result, the V-V distances in decreasing or- 
der are: V2-V3, V1-V2, V1-V4 and V3-V4. To a smaller 
extent, the average V-S distance (taking into account 
only the six nearest-neighbor sulfur ions, respectively) 
also varies. Again, the averaged distance A Jys is largest 
for the V2 ion and smallest for the V4 ion. 

Although the CDW phase is experimentally known^ to 
be paramagnetic insulating, the LDA calculation finds it 
to be metallic (see Fig. [TQ|) . Hence whereas the insuffi- 
ciency of LDA for the metallic regime required a closer 
band-structure study, the failure due to the neglection 
of strong correlations is now obvious. Still the LDA ap- 
proach may deliver relevant information about the kinetic 
part of the hamiltonian and the changes of the electronic 
structure due to the CDW distortion. 

The total LDA DOS for the Im structure is rather 
close to the Cmc2i one. We plotted in Fig. [10] only the 
Wannier bands for the low-energy regime on top of the 
full LDA band structure. The lowest and highest t2g 
bands still have the strongest Aig weight, and the hy- 
bridization between Aig and Egi is significant especially 
at higher energy. However no clear distinction between 
the individual electronic character of the V ions can be 
made on the level of simple projection onto local orbitals. 

A bit more insight is obtained when going to the Wan- 



TABLE V: Onsite terms of the t2g Wannier hamiltonian for 
/m-BaVSs compared to the ones in the metallic phases. The 
values are in meV. 





VI 


V2 


V3 


V4 


Cmc2i 


P63 / mmc 


Aig 


478 


473 


460 


477 


423 


395 


Egl 


431 


431 


436 


432 


210 


200 


Eg2 


448 


444 


443 


448 


236 


200 



TABLE VI: Spread (r^) of the t2g WFs for /m-BaVSs com- 
pared to the ones in the metallic phases. The values are in 
(a.u.2). 





VI 


V2 


V3 


V4 


Cmc2i 


P63 / mmc 


Aig 


22.72 


21.13 


22.18 


23.69 


16.60 


18.62 


Egl 


17.61 


18.04 


18.39 


17.54 


17.55 


17.10 


Eg2 


18.86 


19.31 


17.91 


18.21 


17.53 


17.10 



nier represenation, again using the crystal-field basis as 
the choice of reference. Figure [11] indicates that from 
a low-energy perspective the (V1,V2) and (V3,V4) ions 
have similar characteristics, especially in the occupied 
part of the DOS. The Aig occupation compared to the 
Egl one is larger for the (V3,V4) ions. Because of the 
way of the shifts of the atomic positions in the mono- 
clinic structure it is not so surprising to find that already 
on an LDA level the (V3,V4) and (V1,V2) ions appear to 
form somehow two different classes. Since the (V1,V2) 
ions are more isolated, the larger Egi filling makes sense 
when being the more localized state. In spite of these 
differences it is seen in Tab. |Vl that the LDA crystal- 
field splitting within the t2g manifold is greatly reduced 
in the CDW state. Note also that the spread of the Aig 
WF is significantly enhanced in the LDA description of 



Acivv=+0-10 




FIG. 9: (color online) Single BaVSs chain for the monoclinic 
structure at 40 K. The numbers on the left mark the devi- 
ation of intrachain V-V distance in comparison to the value 
dw=5.37 a.u. in Cmc2i-BaVS3, whereas the numbers on 
the left denote the deviation of the average V-S distance. All 
values in a.u.. The arrows indicate the dominant shifts of the 
ions again with respect to the orthorhombic structure above 
the CDW instability. 
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FIG. 10: (color online) LDA band structure and derived t2g 
Wannier bands (cyan/lightgray) for monoclinic /m-BaVSa- 
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FIG. 11: (color online) LDA DOS based on the derived t2g 
Wannier functions for monoclinic /m-BaVSa- 



/m-BaVSs (see Fig.|Vll). 



D. LDA+DMFT study of the paramagnetic 
insulating regime 

In order to overcome the obvious failure of LDA in 
describing the paramagnetic insulating state of BaVSs 
an LDA+CDMFT approach was employed. Because of 
the revealed different behavior of the inequivalent V ions 
already on the LDA level, the minimal cluster has to 
include all four V ions along the chain. Such a four- 
site cluster leads in the present case to the description 
within an effective twelve-band model on the basis of the 
derived Wannier hamiltonian. The same values for U and 
J as for the metallic regime were used {U =3.5 eV, J=0.7 
eV) , no explicit interatomic Coulomb repulsion term was 
introduced. Due to the large computational effort for the 
QMC impurity solver, the lowest temperature with a still 
reasonable statistics we achieved was T=460 K. 

Figure [12] exhibits the site- and orbital-resolved lo- 



cal spectral function obtained from CDMFT. It is seen 
that within this description the system can be inter- 
preted to be in an insulating state. The partially still 
remaining minor spectral weight at zero energy is due 
to the limitation concerning the handable temperature 
with the QMC solver. Remember that the experimental 
charge gap from precise optics measurements^ equals 
only Ach=0.42 meV. 

Obviously the strong correlations lead to a substantial 
renormalization of the crystal-field splitting, since the Eg 
twofold is now widely separated in energy, shifting the 
Eg2 state towards higher energy and somehow "out of 
the picture". The different tendencies in the character 
of the inequivalent V ions seen in LDA, are now much 
more strongly enhanced. Thus the (V1,V2) ions have 
now nearly exclusively Egi weight, while the (V3,V4) 
ions show some mixed Aig/Egi occupation. This inter- 
esting result is also summarized in the site- and orbital- 
resolved occupations shown in Tab. I VIII Hence the 
charge transfer that was observed in the metallic regime 
due to strong correlations, takes place also in the insu- 
lating state, however now in a site-dependent manner. 
Whereas the (V1,V2) ions loose their Aig occupation 
nearly completely and become orbitally polarized, the 
(V3,V4) ions fall more in the regime of orbital compensa- 
tion, with the V4 ion gaining some extra Egi weight. On 
average, roughly speaking the Egi orbital is the winner of 
the CDW transition, since it replaces the Aig orbital as 
the dominant orbital in the system. Concerning the ques- 
tion of charge order, the data does not provide a strong 
argument for either side. Albeit from the numbers there 
appears to be the slight tendency to put some minor ex- 
tra charge on the (V1,V2) ions, due to the restrictions 
in quantitative accuracy of the formalism this may be 
within the error bars. Remember that experimentally no 
charge order was found^^. 
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FIG. 12: (color online) Local spectral functions from 
LDA+CDMFT for monoclinic /m-BaVSs. The QMC solver 
was used at T=460 K. 
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FIG. 13: (color online) Aig and Egi onsite self-energy J^{icUn) 
for the four inequivalent V ions in monoclinic /m-BaVSs at 
T=460 K. Dark lines correspond to Aig, red/gray lines to 

Egl. 



The LDA+CDMFT method gives access to more quan- 
tities than solely the spectral functions and onsite den- 
sities. To find out more about the nature of the CDW 
transition and the insulating state in BaVSs, investigat- 
ing the behavior of the self-energy Y!,{iuJri) is very in- 
structive. Note that in the present case corresponds 
to a 12x12 matrix, including information not only about 
onsite but also intersite correlations within the four-site 
cluster. 

It is seen in Fig. [13] that the onsite self-energies for 
the various V ions display the expected behavior. While 
TiAig and ^Egi iiave rather different amplitude and 
zero- frequency slope for (V1,V2), they scale similar for 
(V3,V4). The large negative increase of Re ^Egi close 
to zero frequency for (V1,V2) leads to the strong shift 
of the QP spectral weight to lower energies observed in 
Fig. [121 An important qualitative finding is that none of 
the Im E diverges at U0n=0. Hence the opening of the 
gap in BaVSs is due to shifts of the QP states away from 
zero energy. 

The inspection of the nearest-neighbor self-energy 
Ev-v (see Fig. [14]) reveals more details of the CDW 
state. For the (V3,V4) ions Re Sy-v displays a salient 
increasing behavior when approaching zero frequency. 
Such a tendency for the intersite self-energy is a strong 



TABLE VII: Orbital-resolved fillings for /m-BaVSs from 
LDA+CDMFT within the crystal-field Wannier basis. The 
values in the left column correspond to the LDA result, while 
in the right column the values for [7=3.5 eV, J=0.7 eV (with 
T=460 K in the QMC impurity solver) are given. 





VI 


V2 


V3 


V4 


V 


Egl 

Eg2 


0.49 0.12 
0.46 0.89 
0.05 0.03 


0.40 0.11 
0.44 0.85 
0.07 0.07 


0.62 0.47 
0.28 0.46 
0.11 0.03 


0.61 0.34 
0.37 0.62 
0.10 0.02 


0.53 0.26 
0.39 0.70 
0.08 0.03 


sum 


1.00 1.04 


0.91 1.03 


1.01 0.96 


1.08 0.98 





indication for the importance of intersite correlation ef- 
fects leading to interatomic dimer formation^^. How- 
ever within the crystal-field basis (derived from the LDA 
hamiltonian) the occupation on (V3,V4) is of mixed 
(Aig^Egi) character. Diagonalizing the interacting clus- 
ter Green's function thus may lead to an orbital basis 
which corresponds to this dimer symmetry. On the other 
hand, the V1-V2 pair does not show strong intersite cor- 
relations. Thus these dominant Egi occupied ions do not 
tend to form a spin singlet. It follows therefrom that the 
apparent spin gap in the insulating system is not orig- 
inated from the direct correlated coupling of the neigh- 
boring Egl spins on (V1-V2). However, one observes 
for the V2-V3 self-energies a minor tendency for singular 
behavior, especially for the Egi-Aig channel. One may 
conclude from this that the Egi spins on (V1,V2) are 
somewhat effectively pinned by the neighboring dimers 
within the overall tetramerized state. This would ex- 
plain the large drop in the magnetic susceptibility^ and 
excess entropy^ below the MIT. 

It has to be noted that the intersite Sy-v for V4-V1 
is not a true pair self-energy as the others. Since the 
cellular DMFT approach we applied here breaks trans- 
lational invar iance, this special intersite self-energy con- 
necting the surface of the cluster may not be obtained 
accurately. Thus for completeness we plotted in Fig. [14] 
instead Evi-V4- The issue of translational-invariance 
breaking is surely a drawback of the used method. How- 
ever since the V4-V1 pair is expected to be in an interme- 
diate state between VI- V2 and V3-V4, we believe that 
the qualitative result of how the V ions electronically re- 
late to each other remains unchanged. Note that one 
option to restore the periodicity in the present context 
would have been to perform calculations in the chain- 
DMFT^^ framework. However, there the different treat- 
ment of inter- and intrachain hopping might cause other 
problems for BaVSa, since the ID character (as stated) is 
not very strongly indicated from the band hamiltonian. 




1.0 0.5 

« (eV) 



FIG. 14: (color online) Intersite self-energy Sv-v(^^n) for 
the nearest- neighbor V-V pairs in monoclinic /m-BaVSs at 
T=460 K. Note that the shown Ev4-vi is not a true pair 
self-energy here (see text). 
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IV. SUMMARY AND CONCLUSIONS 

The puzzling physics of BaVSs is dominated by the 
competition between the more itinerant Aig state and 
the quasi-locahzed Eg states, which form together the 
t2g manifold of the V(3(i) shell. Some theoretical models 
in the early days tried to rule out one or the other of 
those orbital sectors for playing an essentials role. How- 
ever several recent experimental and theoretical studies 
revealed clearly the importance of the existence of both 
orbital types in order to find a way of understanding the 
complex electronic phases of this system. Thus BaVSs 
appears to be a manifest multiorbital system, and if at 
all some orbital degrees of freedom might just freeze out 
at very low temperatures. 

LDA approaches to the electronic structure of BaVSs 
tend quite naturally to overestimate the itinerant char- 
acter of the system. The Aig associated Wannier orbital 
has an LDA filling of around 70% for the RT hexagonal 
structure, while the inclusion of strong electronic correla- 
tions within the LDA-f-DMFT reduce this filling to about 
50%. This is because a dominant Aig filling becomes just 
too costly in the presence of reasonable mutual Coulomb 
interactions between the electrons. The balanced occupa- 
tion of Aig and Eg seems to be stable when going to the 
orthorhombic structure below Ts, yet the LDA+DMFT 
calculations revealed some tendency towards increasing 
the Eg filling even more. Note however that there is still 
some arbitrariness in what one calls an Aig /Eg orbital^ 
especially in the metallic regime and between different 
phases, and hence some margin in the derived numbers. 
Nonetheless, it was shown that in the paramagnetic in- 
sulating phase with the monoclinic structure the average 
Eg occupation finally reaches about 70%. 

The Eg states show a low QP coherence temperature 
and are expected to be effectively localized for elevated 
temperatures in the metal. Though the LDA approach 
(even when extended by linear expansions of the DMFT 
self-energ}^^) shows Eg FS sheets, it is very likely that 
the Eg QPs do not participate in the true Fermi surface in 
a well-defined manner. Hence due to the rather large dif- 
ference in the coherence temperatures, one may describe 
BaVSs to be in an effective orbital-selective insulating 
regime. As there is substantial hybridization between 
Aig and Egi below Ts, the low metallicity of BaVSs, 
and especially the bad- metal regime below 150 K, may re- 
sult from the scattering processes for the (quasi-) itinerant 
electrons. The origin of the hexagonal- to-orthorhombic 
transition appears to be closely related to the Aig-Egi 
hybridization. It was shown that the Egi Wannier orbital 
connects to two S2 ions, while the Eg2 one hybridizes with 
only one apical SI ion. It is likely that this imbalance 
favors the susceptibility for growing symmetry-breaking 
^ig-Egi hopping. Thereby the large DOS close to the 
Fermi level is partly reduced when lifiting the Eg degen- 
eracy, and this driving force leads via the final zigzag in- 
trachain distortion to a new energetic minimum. Albeit 
many details of the MIT are still open, as discussed in 



Ref. [13 the correlation-induced fiattened Aig FS sheets 
away from the high-symmetry directions in the BZ are 
good candidates for a meaningful matching with exper- 
imental findingSii^. Further experimental studies of the 
low-energy regime close to the MIT are needed to reveal 
more details. 

We revealed with LDA+CDMFT calculations that the 
insulating CDW state does not stay behind in terms 
of complexity of the electronic structure in comparison 
with the metallic regime. The tetramerization of the 
V ions, structurally an effective trimerization, leads to 
quite different behavior. While the (V3,V4) pair ap- 
parently forms a correlated dimer with mixed Aig /Egi 
occupation, the (V1,V2) ions are strongly orbitally po- 
larized with major Egi occupation and negligible cou- 
pling. Note however that this picture is of course basis 
dependent, and we worked always in the Wannier basis 
derived from the LDA hamiltonian. Hence a new or- 
bital basis may be found in the interacting regime where 
e.g. the occupation for (V3,V4) has also polarized char- 
acter. Concerning the "free" spins on (VI, V2), note that 
they are still coupled to the neighboring dimers, whereby 
their degrees of freedom are substantially reduced. The 
latter effect may serve as an explanation of the quenched 
local moments observed below Tmit • Fagot et al}^ pro- 
posed from anomalous x-ray measurements a dominant 
Egi occupation on VI and an additional dominant oc- 
cupation of Aig on V3, as well as no definite preferen- 
tial occupation on (V2,V4). This would describe an or- 
bital order modulated with 2c along the chain. Our pic- 
ture differs by the fact that we do not find a dominant 
Aig occupation, but rather two ions, i.e., (V1,V2), with 
dominant Egi weight. Although from a local structural 
point of view the former proposition appears meaning- 
ful, the intersite correlation effects appear efficient in sin- 
gling out dimer and isolated behavior. Note that dimer 
formation and nearest-neighbor spins may also be ener- 
getic favorable, since the former brings in some residual 
hopping whereas the latter yields energy from spin ex- 
change. Recently, Fazekas et al.^^ proposed a minimal 
one-dimensional model, expecting thereby also nearest- 
neighbor Eg spins. However the orbital degeneracy of the 
Eg multiplet, ascribed to spin-orbit coupling, was kept in 
that model. 

The aspect of the spin degree of freedom in BaVSs 
was so far only scarcely discussed. Albeit it is expected 
that the former is an important ingredient in the under- 
standing of the physics, its role in the different phases 
is still quite open. Nakamura and coworkersi^ reported 
anomalies in the temperature derivative of the magnetic 
susceptibility not only at Tmit but also at the onset of 
the bad- metal regime (^160 K). Hence the detailed role 
of the local Eg spins in the metal still poses important 
questions. Our study of the insulator renders it difficult 
to make precise statements about possible long-range or- 
der effects for the spin and orbit als, thus leaving ques- 
tions concerning the classification of the electronic phase 
open. Furthermore, the mechanisms in conjunction with 
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the famous magnetic transition at Tx need to be ad- 
dressed. LDA+(C)DMFT computations may still help 
in delivering some further information associated with 
those problems. For instance, the calculation of mag- 
netic susceptibilities in the insulating regime, the study of 
exchange interactions, as well as the investigation of the 
influence of the interchain coupling especially on the spin 
arrangement. More sophisticated model studies for mul- 
tiorbital chains, perhaps with therefore well-adapted the- 
oretical tools like density matrix renormalization group 
(DMRG)^, may also deliver important new insights. 
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